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Stability Analysis of Slotted Aloha with 
Opportunistic RF Energy Harvesting 

Abdelrahman M.Ibrahim, Ozgur Ercetin, and Tamer ElBatt 


Abstract 

Energy harvesting (EH) is a promising technology for realizing energy efficient wireless networks. In 
this paper, we utilize the ambient RE energy, particularly interference from neighboring transmissions, to 
replenish the batteries of the EH enabled nodes. However, RE energy harvesting imposes new challenges 
into the analysis of wireless networks. Our objective in this work is to investigate the performance of a 
slotted Aloha random access wireless network consisting of two types of nodes, namely Type I which 
has unlimited energy supply and Type II which is solely powered by an RE energy harvesting circuit. The 
transmissions of a Type I node are recycled by a Type II node to replenish its battery. We characterize 
an inner bound on the stable throughput region under half-duplex and full-duplex energy harvesting 
paradigms as well as for the finite capacity battery case. We present numerical results that validate our 
analytical results, and demonstrate their utility for the analysis of the exact system. 

Index Terms —Wireless networks, slotted Aloha, opportunistic energy harvesting, interacting queues. 


I. Introduction 

One of the prominent challenges in the field of communication networks today is the design of 
energy efficient systems. In traditional networks, wireless nodes are powered by limited capacity 
batteries which should be regularly charged or replaced. Energy harvesting has been recognized as 
a promising solution to replenish batteries without using any physical connections for charging. 
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Nodes may harvest energy through solar eells, piezoeleetrie deviees, RF signals, ete. In this 
paper, we foeus on RF energy harvesting. Reeent studies present experimental measurements 
for the amount of RF energy that ean be harvested from various RF energy sourees. Two main 
faetors affeet the amount of RF energy that ean be harvested, namely, the frequeney of the RF 
signal and the distanee between the “interfered’ and the harvesting node, e.g., see Table I in [[T]|. 

Reeently, an information-theoretie study of the eapaeity of an additive white Gaussian noise 
(AWGN) ehannel with stoehastie energy harvesting at the transmitter has shown that it is equal to 
the eapaeity of an AWGN ehannel under an average power eonstraint |l2l. This, in turn, motivated 
the investigation of optimal transmission polieies Q for single user Hj-jbll and multi-user fTlI- 
ll9l energy harvesting networks. The optimal poliey that minimizes the transmission eompletion 
time was studied in dU. In [l5]|, the authors studied the problem of maximizing the short-term 
throughput and have shown that it is elosely related to the transmission eompletion time problem 
ll4| . The authors in |l6l studied the optimal transmission polieies for energy harvesting networks 
under fading ehannels. Moreover, [|7]|-[|9l| extends the analysis to broadeast, multiple aeeess, 
and interferenee ehannels, respeetively. The authors in ifTOll introdueed the eoneept of energy 
eooperation where a user ean transfer portion of its energy, over a separate ehannel, to assist 
other users. 

Signifieant researeh has also been eondueted on RF energy harvesting. In [fTTll . the author 
diseusses the fundamental trade-offs between transmitting energy and transmitting information 
over a single noisy link. The author derives the eapaeity-energy funetions for several ehannels. 
The authors in ifT^ . extend the point-to-point results of ifTTI to multiple aeeess and multi¬ 
hop ehannels. Reeently, several teehniques were proposed for designing RF energy harvesting 
networks (RF-EHNs), e.g., ifBl . The RF energy harvesting proeess ean be elassified as follows: 
i) Wireless energy transfer, where the transmitted signals by the RF souree are dedieated to energy 
transfer, ii) Simultaneous wireless information and power transfer, where the transmitted RF 
signal is utilized for both information deeoding and RF energy harvesting, and Hi) Opportunistic 
energy harvesting, where the ambient RF signals, eonsidered as interferenee for data transmission, 
are utilized for RF energy harvesting. The reeeiver arehiteeture may also vary as follows llT3]l . 
HU: a) Co-located receiver architecture, where the radio reeeiver and the harvesting eireuit 
use the same antenna for both deeoding the data and energy harvesting, and b) Separated 
receiver architecture, where the radio reeeiver and the harvesting eireuit are separated, and 
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each is equipped with its own antenna and RF front-end circuitry. 

In HU, the authors discuss the practical limitations of implementing a simultaneous wireless 
information and power transfer (SWIFT) system. A major issue is that energy harvesting circuits 
are not able to simultaneously decode information and harvest energy. Hence, the authors in lfT4]| 
proposed and analyzed two modes of operation for the co-located receiver architecture, that is, 
time switching and power splitting. Furthermore, the RF energy harvesting transceivers may 
also be classified as: a) Half-duplex energy harvester, where a co-located RF energy harvesting 
transceiver can either transmit data or harvest RF signals at a given instant of time, and b) Full- 
duplex energy harvester, where a node is equipped with two independent antennas and can 
transmit data and harvest RF signals, simultaneously. In this paper, we investigate opportunistic 
RF energy harvesting under the half-duplex and full-duplex modes of operation. 

The cornerstone of random access medium access control (MAC) protocols is the Aloha 
protocol ifTSll . which is widely studied in multiple access communication systems because of its 
simplicity. The applications of Aloha-based protocols range from traditional satellite networks 
[fT^ to radio frequency identification (RFID) systems [fTTll and the emerging Machine-to Machine 
(M2M) communications [fT^ . It is also considered as a benchmark for evaluating the performance 
of more sophisticated MAC protocols. Based on the Aloha protocol, nodes contend for the 
shared wireless medium and cause interference to each other. Hence, the service rate of a node 
depends on the backlog of other nodes, i.e., the nodes’ queues become interacting as originally 
characterized in |[T9l . Tsybakov and Mikhailov |[20]| were the first to analyze the stability of 
a slotted Aloha system with finite number of users. Rao and Ephremides [|2TI characterized 
the sufficient and necessary conditions for queue stability of the two user case, using the so- 
called stochastic dominance technique. In addition, they established conditions for the stability 
of the symmetric multi-user case. Other works followed and studied the stability of slotted Aloha 
with more than two users Il22ll - ll25]l . The authors in [l26l extended the stability analysis under the 
collision model to a symmetric multi-packet reception (MPR) model, which was later generalized 
to the asymmetric MPR model in ^TT\ . 

Perhaps the closest to our work are [|28]| . |[29]l which characterize the stability region of a slotted 
Aloha system with energy harvesting capabilities, under the multi-packet reception model. The 
authors considered a system where the nodes harvest energy from the environment at a fixed 
rate and, thus, the energy harvesting process is modeled as a Bernoulli process. 
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In this work, we analyze the stability of a slotted Aloha random aeeess wireless network 
eonsisting of nodes with and without RF energy harvesting eapability. Speeifieally, we eonsider 
a wireless network eonsisting of two nodes, namely a node of Type I whieh has unlimited energy 
supply and a node of Type II which is powered by an RF energy harvesting circuit. The RF 
transmissions of the Type I node are harvested by the Type II node to replenish its battery. 
Our contribution in this paper is multi-fold. First, we outline the difficulties in analyzing the 
stability of the exact RF energy harvesting Aloha system So and for mathematical tractability we 
introduce an equivalent system So- Second, we generalize the stochastic dominance technique 
for analyzing RF EH-networks. Third, we characterize an inner bound on the stable throughput 
region of the system So under the half-duplex and full-duplex energy harvesting paradigms. 
Also, we derive the closure of the inner bound over all transmission probability vectors. Fourth, 
we investigate the impact of finite capacity batteries on the stable throughput region. Finally, 
we validate our analytical findings with simulations and conjecture that the inner bound of the 
system Sq is also an inner bound for the exact system Sq- 

The rest of this paper is organized as follows. In Section HU we present the system model 
and the assumptions underlying our analysis. In Section Hill we describe the energy harvesting 
models for the systems Sq and So- Our main results are presented in Section HY] and proved in 
Section |Vl In Section |Vll we investigate the impact of finite capacity batteries and full-duplex 
energy harvesting on the stability region of our system. We corroborate our analytical findings 
by simulating the systems Sq and So in Section IVIII Finally, we draw our conclusions and point 
out directions for future research in Section IVllli 

II. System Model 

We consider a wireless network consisting of two source nodes and a common destination, as 
shown in Fig. [T1 We consider a slotted Aloha multiple access channel [flSl . where time is slotted 
and the slot duration is equal to one packet transmission time. We assume two types of nodes in 
our system: Type I node has a data queue, Qi, and unlimited energy supply, while Type II node 
has a data queue, Q 2 , and a battery queue, B, as shown in Fig. [T] Moreover, packets arrive to 
the data queues, Qi and Q 2 , according to independent Bernoulli processes with rates Ai and A 2 , 
respectively. The transmission probabilities of Type I and II nodes are qi and q 2 , respectively. 

We assume perfect data channels, i.e., the destination successfully decodes a data packet. 
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Type I 



Fig. 1: System model 


if only one node transmits. If two nodes transmit simultaneously, a eollision oeeurs and both 
paekets are lost and have to be retransmitted in future slots. At the end of eaeh time slot, the 
destination sends an immediate aeknowledgment (ACK) via an error-free feedbaek ehannel. 

Data paekets of Type I and II nodes are stored in queues Qi and Q 2 , respeetively. The evolution 
of the queue lengths is given by iflSll 

Q*+i=max{Q*-r/,0} + X^ * = 1,2. (1) 

where Xj G {0,1} is the arrival proeess for data paekets and G {0,1} is the departure proeess 
independent of Qi status, i.e., F/ = 1 even if the data queue is empty ETll . X* is a Bernoulli 
proeess with rate A* and E[y/] = /x*. 

We assume that Type II nodes operate under a half-duplex energy harvesting mode, i.e., they 
either harvest or transmit but not both simultaneousljQ. Henee, the harvesting opportunities are 
in those slots when a Type II node is idle while a Type I node is transmitting. The ehannel 
between the two souree nodes is a bloek fading ehannel, where the fading eoeffieient remains 
eonstant within a single time slot and ehanges independently from a slot to another. For Rayleigh 
fading, the instantaneous ehannel power gain ht at time slot t is exponentially distributed, i.e., 
ht ~ Exp{l). Let Pj be the transmission power of node j, and / be the distanee between the 
two souree nodes. We eonsider the non-singulaj^ pathloss model with (1 -f where a is 

*We extend the analysis to the case of full-duplex RF EH in Section Ivll 

^In order to harvest significant amount of RF energy, I is typically small. Hence, we use a non-singular (bounded) pathloss 
model instead of a singular (unbounded) pathloss l~°‘ model, because the singular pathloss model is not correct for small values 
of I due to singularity at 0, 1301 . 
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the pathless exponent. We use power and energy interehangeably throughout the paper, since 
we assume unit time slots. In the next section, we develop a discrete-time stochastic process to 
model RF energy harvesting. 


III. Energy Harvesting Models 


Prior work largely models the energy harvesting process as an independent and identically 
distributed (iid) Bernoulli process with a constant rate ll29ll^ . In the following, we first study the 
RF energy harvesting process of the exact system Sq and show that the inter-arrival time of the 
energy arrivals has a general distribution. Second, we propose an equivalent system Sq where 
the inter-arrival time of the energy arrivals is geometrically distributed with the same mean as 
5o. 


A. RF Energy Harvesting Model of the exact system So 

The received power at a Type II node from the transmissions of Type I node at time slot t is 
Pnit) = where Pi is the transmission power of Type I node and 77 is the RF 

harvesting efficiency [fTSll . Recent studies demonstrated that 77 typically ranges from 0.5 to 0.7, 
where its value depends on the efficiency of the harvesting antenna, impedance matching circuit 
and the voltage multiplier [1311 . In order to develop the analytical model underlying this paper, 
we approximate the continuous energy arrival process in quantas of size 7 joules. 

Typically, we need to harvest RF energy from multiple transmissions of Type I node in order 
to accumulate 7 joules. Conceptually, accumulating 7 joules is equivalent to having a single 
energy packet arrival to the battery. We model the battery of a Type II node as a queue with 
unit energy packet arrivals from the harvesting process, and unit energy packet departures when 
Type II node transmits. There is an energy packet arrival to the battery queue at the end of the 
time slot in which the accumulated energy exceeds 7, see Fig. |2l Let Z be the number of Type 
I transmissions needed to harvest one energy unit. 

Lemma 1. For a persistent (gi = l) and saturated Type I node, the probability mass function 
(PMF) of Z, when the channel between Type I and II nodes is modeled as a Rayleigh fading 


^Typically, the energy harvesting process is not iid, because to harvest one energy unit, energy is accumulated over multiple 
slots. 
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Fig. 2: A unit of energy is harvested every time the aeeumulated energy in the battery exeeeds 

7 - 


ehannel with parameter ht, is given by 

-e nk-l /I I ia\ 

P[Z = fc] = ^^, = where0 = ^i^ (2) 

Proof: Note that if the aeeumulated reeeived energy over k slots is greater than or equal to, 
7 while the aeeumulated reeeived energy up to slot fc — 1 is less than 7, then we need k slots 
to harvest one energy unit. 
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The distribution of sum of independent exponential random variables is an Erlang Distribution. 
Henee, ~ Erlang{k,l), where k is the shape parameter and ht ~ Exp{l). From, the 
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cumulative distribution function (CDF) of the Erlang distribution, we know that 


P 
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1 _ e-® 


j=0 




(8) 


Hence, (b) is obtained by substituting ([8]) in ([4]). ■ 

From (|5]), we notice that the distribution of the inter-arrival times of the energy harvesting 
process is a shifted Poisson distribution, i.e., Z = V + 1, where E is a Poisson random variable 
with mean 6. The expected inter-arrival time of the harvesting process is given by E[Z] = 1-1-6'. In 
the general case where Type I node is unsaturated and transmits with probability qi, characterizing 
the PMF of Z is challenging because the queue evolution process is not an iid process. For 
instance, if Q\ = 2, we know for sure that > 1. 


B. RF Energy Harvesting Model of the equivalent system So 

For the mathematical tractability of the results derived in later sections, we consider a system 
Sg, where the RF energy harvesting process is an iid Bernoulli process. Fet Ph\{i} be the mean 
of the iid Bernoulli process of unit energy packet arrivals, where it can be interpreted as the 
probability of success in harvesting one energy unit given that Type I node is transmitting. Based 
on the fact that the inter-arrival time of a Bernoulli process is geometrically distributed, the mean 
inter-arrival time is l/p/i|{i}. Hence, the relationship between the exact harvesting process in Sq 
and the equivalent Bernoulli process in Sq is given by 


PHUD = 


(9) 


In general, an iid Bernoulli process has a rate Ph\M^ where Ph\M is the probability of harvesting 
one energy unit given a set of nodes M. are transmitting. Under half-duplex energy harvesting. 
Type II node only harvests from the transmissions of Type I node, when the Type II node is not 
transmitting, i.e., the probability of harvesting one energy unit given Type II node is transmitting 
Ph\{ 2 } = 0, and the probability of harvesting one energy unit given both nodes are transmitting 
Ph\{i, 2 } = 0. For convenience, we denote Ph\{i} by ph- 


C. Analyzing the Battery queue in Sq 

We assume that a Type II node opportunistically harvests RF energy packets from the trans¬ 
missions of Type I node. Also, transmitting a single data packet costs one energy unit. Fet 
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denotes the energy harvesting process modeled as a Bernoulli process. Assuming half-duplex 
harvestingjj the average harvesting rate is the difference between the fraction of time slots in 
which the Type I node is transmitting and the fraction of time slots in which both nodes are 
transmitting, i.e., 

E[H^] = qiph P[Qi>0] - qiq2Ph P[Qi > 0, 5 > 0, Q 2 >0]. (10) 


The battery queue evolves as 


Qt+l ^ 


( 11 ) 


where /i^ G {0,1} represents the energy consumed in the transmission of a data packet at time 
t. Under backlogged data queues Qi and Q 2 , the average rate of harvesting becomes 

E[H^\Qi>Q,Q 2>Q] = qiPh{l - q2E[B>Q\Qi>Q,Q2>Q]Y ( 12 ) 


Now, the energy harvesting rate is only dependent on the battery queue status. Thus, if the 
battery queue is empty, the energy harvesting rate is qiPh, otherwise it is reduced to qiPh{^ — <12) 
because of the half-duplex operation. Hence, the battery queue forms a decoupled discrete-time 
Markov chain, as shown in Fig. |3l By analyzing the Markov chain we find the probability that 
the battery is non-empty, as given by the following lemma. 

Lemma 2. For a half-duplex RF energy harvesting node with infinite capacity battery, the 
probability that the battery is non-empty, given that the data queues Qi and Q 2 are backlogged, 
is given by 


F[B > 0|Qi > 0, Q 2 > 0] = min 


qiPh 

^ 2(1 + qiPh) 



(13) 


Proof: Let tt = [tto, tti, ...] be the steady-state distribution of the Markov chain shown in Fig. 
Ill Applying the detailed balance equations, we obtain j (1 — ?2)^~^7ro, /c = 1, 2, • • • . 

Therefore, by substituting in the normalization condition tt* = 1, we get the utilization factor 
p = 1 - TTo = Hence, F[B > 0|Qi > 0, > 0] = min {p, 1} . ■ 

IV. Main Results 


In this section, we present our main results pertaining to the stable throughput region of the 
opportunistic RF energy harvesting slotted Aloha network So- We adopt the notion of stability 


''in the sequel, we will extend the model to full-duplex as well. 
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Fig. 3: Markov chain model of the battery queue given that the data queues Qi and Q 2 are 
baeklogged. Note that 5' = l—qiph— <?2 and 5 = l — qiph{^ — q' 2 )— <? 2 - 


proposed in [|^ . where the stability of a queue is determined by the existenee of a proper 
limiting distribution. A queue is said to be stable if 


lim P[Q* < x] = F{x) and lim F{x) = 1. (14) 

t —>-00 iT —>-00 


The stability of a queue is equivalent to the reeurrence of the Markov chain modeling the 
queue length. Loynes’ Theorem [133 states that if the arrival and service proeesses of a queue 
are strietly jointly stationary and the average arrival rate is less than the average serviee rate, 
then the queue is stable. Also, if the average arrival rate is greater than the average service rate, 
then the queue is unstable and the queue size approaehes infinity almost surely. The stable 
throughput region of a system is defined as the set of arrival rate veetors, (Ai, A 2 ) for our system, 
for which all data queues in the system are stable. 

Next, we establish suffieient conditions on the stability of the opportunistie RF energy har¬ 
vesting Aloha Sg- Assuming half-duplex energy harvesting and unlimited battery capaeity, the 
stability region is characterized by the following theorem. 

Theorem 1. An inner bound on the stable throughput region of the opportunistic RF energy 
harvesting slotted Aloha Sg is the triangle OBD, shown in Fig. HI Assuming half-duplex energy 
harvesting and unlimited battery size, the region is charaeterized by 


pinner 


— \ (Ai, A 2 ) Ai < gi ( 1 


Ao 


1 - gi 


, A 2 < 




gi (1 — min 


giPh 

(i+qiph) 


,?2 


(15) 


Proof: The proof is established in the Seetions IV-AI to IV-El 


Theorem 2. The elosure of the inner bound over all transmission probability veetors 
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q = {qi,q 2 ) is characterized by 



U S.(9) 


(Ai,A2) A 2 < — ^1 — Ai + A2+\/(1+Ai —A2)^— 4Ai^ 


Proof: The proof is established in Seetion IV-F[ 


( 16 ) 


V. Stability Analysis 

For the majority of prior work on stability analysis of interaeting queues, the serviee rate of a 
typieal node deereases with respeet to the transmissions of other nodes in the system. Perhaps, 
the most basie example is the eonventional slotted Aloha system [lUll . where inereasing the 
serviee rate of an arbitrary node eomes at the expense of deereasing the serviee rate of other 
nodes. For our purposes, we refer to sueh systems without energy limitations as “interferenee- 
limited” systems. 

On the other hand, in our RF energy harvesting system, transmissions from interfering nodes 
give rise to two opposing effeets on Type II (RF energy harvesting) nodes. Similar to elassie 
interferenee-limited systems, the interfering nodes ereate eollisions and, thus, deerease the ser¬ 
viee rate of RF energy harvesting nodes. Meanwhile, transmissions from interfering nodes are 
exploited by RF energy harvesting nodes to opportunistieally replenish their batteries. Therefore, 
from the perspeetive of an RF energy harvesting node, a fundamental trade-off prevails between 
the inereased number of energy harvesting opportunities and the inereased eollision rate, whieh 
are both eaused by interferenee. As will be shown formally, this fundamental trade-off splits the 
stable throughput region for RF energy harvesting slotted Aloha networks into two sub-regions, 
a sub-region where interferenee is advantageous for the RF energy harvesting node and another 
sub-region where it is not. These two sub-regions map direetly to two modes of operation for 
our system and are eharaeterized as follows: 

1) Energy-limited mode: is the sub-region of the stable throughput region in whieh the trans¬ 
missions of interfering nodes enhance the throughput of the RF energy harvesting node, 
i.e, the throughput enhaneement due to the inereased harvesting opportunities outweights 
the degradation due to eollisions ereated by the interfering nodes. 

2) Interference-limited mode: is the sub-region of the stable throughput region in whieh the 
transmissions of interfering nodes do not increase the throughput of the RF energy har- 
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Fig. 4: The energy limited and interference limited sub-regions of the stable throughput region 
of Sn are characterized by the triangles BCO and BCD, respectively. 


vesting node, i.e, the throughput degradation due to collisions equals or outweighs the 
throughput increase due to the increased energy harvesting opportunities. 

The two parts of the stable throughput region for our system are shown in Fig. IH where the 
energy-limited region is enclosed by the triangle BCO and the interference-limited region is 
enclosed by the triangle BCD. 

In order to characterize the inner bound on the stability region of Sq in Theorem 1, we 
introduce a deprived system Sn where Type II node transmits only in the time slots in which 
Qi is non-empty. We derive the stability region of Sd, by going through the following three 
steps discussed next. First, we characterize the average service rates of the two interacting data 
queues. Second, we generalize the Stochastic dominance technique proposed in [1211 to capture 
our system dynamics and two-mode operation. Third, we derive the stability conditions of 
using the generalized stochastic dominance approach. Finally, we prove in Section IV-EI that the 
stability region of is an inner bound on the stability region of Sq- 

A. Service Rates of the Interacting Queues in Sq 

The average service rates of the data queues, Qi and Q 2 , in Sg are given by 

/ii = gi(l-g2 P[ 5 > 0 ,g 2 > 0 |gi> 0 ]), (17) 

fi2 = <12 ^[B >0, Qi = 0\Q2>0]q2 — Qi) P[i?>0, gi >0|g2>0]. (18) 
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where the serviee rate of the Type I node, Hi, is the fraction of time in which Type I node 
decides to transmit, excluding the fraction of time in which Type II node is also transmitting. A 
Type II node transmits, if it is active, i.e., i?>0 and Q2>0, and decides to transmit. Similarly, 
the service rate of the Type II node, yU2, is the fraction of time in which Type II node has a 
non-empty battery and decides to transmit, excluding the fraction of time in which both nodes 
are transmitting. Note that the queue evolution equation in O implies that P[(5i > 0] = 1 in 
(fTTl) and F[Q 2 > 0] = 1 in (fTSl) . 

In our system, we have three interacting queues, namely Qi, Q 2 and B. The analysis of three 
interacting queues is prohibitive and, hence, calculating the probability P[i? >0, Qi = 0|(52 >0]. 
Therefore, we consider a system Sd where ¥[B > 0,Qi = 0\Q2 > 0] = 0, i.e., we consider a 
lower service rate for the Type II node, since it transmits only in the time slots in which Qi is 
non-empty. The relationship between the two systems is discussed in Section IV-EI 

B. Service Rates of the Interacting Queues in So 

In order to analyze the interaction between Qi and Q 2 in Sd, we decouple the battery queue, 
B, by substituting the probability of the battery queue being non-empty with the conditional 
probability given by Lemma 2. Hence, the average service rates of the data queues Qi and Q 2 


are given as = gi (l - ga P[5 >0|Qi > 0, Qa > 0] P[Q 2 >0|Qi > 0]) (19) 

= gi(^l - min P[Q 2 > 0 |Qi> 0 ]j. (20) 

fi2 = q2{l-qi)^B>0\Qi>0,Q2>0] P[gi>0|g2>0] (21) 

= (1 -gi)min|^^-^i^^,g2} P[gi > 0|g2 >0]. (22) 


The probability that Qi is non-empty given that Qj is saturated (always backlogged) is given by 
P[gi > 0|gj > 0] = ^, z = 1, 2, z 7^ j, where /z| is the service rate of Qi given that both data 
queues are saturated. We derive z = 1,2 in Section |WDl 

From (I20I) . we note that /ii decreases with increasing P[g2 > o|gi > 0]. Recall that for 
the interference-limited region, increasing the service rate of one node comes at the expense of 
decreasing the service rate of other nodes. The system is interference-limited from the perspective 
of Type I node, since increasing A 2 comes at the expense of decreasing /zi. Also, from (l22l) . we 
observe that ^2 is increasing in P[gi > 0|g2 > 0]. Hence, increasing Ai increases /Z2 until both 
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data queues are saturated. Thus, the system is energy-limited from the perspeetive of the Type 
II node until the data queues beeome saturated. 

In Fig.m we depiet the stable throughput region R/j of the system Sjy. The boundary between 
the energy-limited part and the interferenee limited part, from the perspeetive of the Type II node 
is Ai = Ai, where Ai is the arrival rate for Qi at whieh both data queues, Qi and Q 2 , beeome 
saturated. Aeeordingly, the energy-limited sub-region is eharaeterized by 



(23) 


and the interferenee limited sub-region is eharaeterized by 



(24) 


C. The Generalized Stochastic Dominance Approach 

In this seetion, we rely on stoehastie dominanee arguments iflTI . whieh are instrumental in 
establishing the stable throughput region of So- However, the eonventional stoehastie dominanee 
approaeh should be modified for our system, sinee the transmission of dummy paekets by Type 
I node inereases the harvesting opportunities for Type II node. Thus, in order to eonstruet 
a hypothetieal “dominant” system in whieh the queue lengths are never smaller than their 
eounterparts in the system So, the hypothetieal system proposed in ETI is modified to eapture 
the two-mode operation inherent to our RF EH system. 

Reeall, from (l22l) . that the serviee rate of Q 2 inereases with Ai. Thus, it is straightforward to 
show that, using elassie stoehastie dominanee arguments, saturating Qi inereases the serviee rate 
of Q 2 - Henee, the queue length in this hypothetieal system (partieularly Q 2 ) no longer dominates 
(i.e. eould be smaller) its eounterpart in the system So and, thus, the elassie stoehastie dominanee 
argument fails. For example, if we eonsider the ease where Ai = 0 and A2 > 0, we observe that 
A2 < Q'2 belongs to the stable throughput region, whieh eontradiets (l22l) . where for Ai = 0, we 
get A2 = 0. 

Henee, we define the hypothetieal system for our RF EH system to be eonstrueted as follows: 

• Arrivals at data queues Qi and Q 2 oeeur at the same instants as in the system So- 

• Transmission deeisions, determined by independent eoin tosses, are identieal to those in the 


system So- 
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• In the energy-limited region, Type I node does not transmit dummy packets. Thus, the 
energy arrivals to the battery queue of the harvesting node occurs exactly at the same 
instants as in the system Sd- 

• In the interference-limited region, if the data queue is empty and the node decides to 
transmit, a dummy packet is transmitted, if the node has sufficient energy to transmit. 

From the construction of the new hypothetical system proposed above, it can be noticed that 
it behaves like the system Sd in the energy-limited region and dominates the system Sd in the 
interference-limited region. Thus, this hypothetical system dominates our system Sd because 
the transmissions of dummy packets collide with the transmission of the other node. Also, the 
transmissions of dummy packets consume energy without contributing to the throughput of Type 
II node. 


D. Establishing the Stability Conditions of Sd 


In order to derive the stability conditions of the system Sd, we construct two dominant 
systems, where in the first dominant system Type II node is backlogged and in the second 
dominant system. Type I node is backlogged only in the interference-limited region. 

1) First dominant system: In this hypothetical system, we consider the case where the Type 
II node continues transmitting dummy packets whenever its data queue, Q 2 , is empty given that 
its battery is non-empty. Since the system is interference-limited from the perspective of Type I 
node, our dominant system is identical to the one proposed in ETI . Hence, the saturated service 
rate of Qi is given by 

Also, by substituting P[(5i > OIQ 2 > 0] = Ai//i* in (l22l) . the service rate of Q 2 becomes 

d2 = -7- -f -YV- 


Therefore, the stable throughput region of the first dominant system Si is given by 


!. = I (Ai. A,) I A. < „ (1 - mm {fc}). A, < 


(l-ft)mm{(lg|j.<i2|A 


(27) 
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Also, since the system beeomes interferenee-limited from the perspective of node II, when 
the data queues, Qi and Q 2 , are baeklogged, Ai is given by 

2) Second dominant system: In this hypothetical system, Qi is baeklogged only in the interferenee- 
limited region from the perspective of Type II node. In the interferenee-limited part of M 2 > the 
saturated serviee rate of Q 2 is given by 

Similarly, by substituting F[Q 2 > 0|Qi > 0] = X 2 /dl we obtain 

= gi (1 - ■ (30) 

Therefore, the stable throughput region of the interferenee-limited part of the second dominant 
system is given by 

K 2 |Ai>Ai = {(Ai,A 2 ) I Ai < Ai < gi(^l - y^)A 2 < (1 - gi) min ^ 2 } }• 

(31) 

The stable throughput regions of the dominant systems Mi and are shown in Fig. HI 

3) Stability region of the system Sd-' In the following lemma we derive the relationship 
between the stable throughput region of the dominant systems Mi and ®2|Ai and the original 
system Rjj. 

Lemma 3. The stable throughput region M^i of the system is given by the union of the 
stable throughput region of the first dominant system and the interference-limited part of the 
second dominant system, i.e., M^^ = Mi LJ ^2|Ai>Ar 

Proof: The stable throughput region of the original system is the union of the two dominant 
systems, based on II2TI . i.e., M^i = M 1 UM 2 . From the construetion of the seeond dominant system, 
we know that the energy-limited region is identical to the original system, i.e., ®2|Ai<Ai ~ 
®z)|Ai<Ai- Henee, we have M 2 = ^d\x,<m ^ ^ 2 |Ai>Ai and Mz, = Mi U U 

Now, assume that the rate pair (xi,X 2 ) G Henee, xi < Ai and X 2 < /U. 2 . Sinee we 
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achieve the maximum service rate for Type II node by backlogging Q 2 , from (1271) we obtain 




(32) 


Therefore, the rate pair {xi,X 2 ) G Mi and ^D|Ai<Ai - 

Proposition 1. The stability region of the system is given by 




— \ (-^ 1 , ^2) 


Ai < gi 1 





(33) 


Proof: For the purpose of the proof, we also define a conventional Aloha system with nodes 

with unlimited energy supplies, i.e., iSAioha- The outline of the proof is as follows 

(i) The generalized stoehastic dominance approaeh proves the neeessity of the stability eon- 
ditions on Qi and Q 2 in (1331) . Meanwhile, it only proves the suffieieney of the stability 
conditions for Q 2 - 

(ii) The stability eondition on Qi in iSAioha is suffieient for stability of Qi in So- 

(iii) The stability eondition on Qi is the same in both systems iSAioha and So- 

(iv) From (ii) and (iii), we establish the suffieieney of the stability eondition on Qi for So- 

The detailed proof is as follows 

(i) Reeall that for systems with unlimited energy sueh as iSAioha, the suffieient and neeessary 
stability eonditions are given by the union of the stability regions of the two hypothetieal 
systems in ETI . On the other hand, for a system with batteries the transmission of dummy 
paekets in the hypothetieal system wastes the energy of the nodes whieh limits the data 
transmissions in future slots. For instanee, in a hypothetieal system where the first node 
is backlogged, there may exist instants at which the first node is unable to transmit due 
to energy outage, while it is able to transmit in the original system. Henee, the second 
node in the hypothetieal system may have a higher sueeess rate eompared to the original 
system, and the sample path dominance is violated. Consequently, the union of the stability 
eonditions of two hypothetieal systems is only a neeessary eondition for the stability of the 
original system ll29ll . In our paper, we arrived to (l33l) by applying the generalized stoehastie 
dominanee approaeh, in whieh the first dominant system is eonstrueted sueh that Type II 
node eontinues transmitting dummy paekets whenever its data queue, Q 2 , is empty. The 
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transmission of dummy packets from Type II node in this system wastes energy. Therefore, 
Type I node in this dominant system may have a better sueeess rate eompared to that of 
the system Sn- Therefore, the generalized stoehastie dominanee approaeh only proves the 
neeessity of the stability eondition on Qi in (|3^ . On the other hand, the stability eondition 
on Q 2 is suffieient and neeessary, sinee Type I node has unlimited energy and the previous 
argument only applies to nodes with batteries. 

(ii) In iSAioha the two eontending nodes have unlimited energy, while in the transmissions of 
the Type II node are eonstrained by the energy in the battery. Henee, the seeond eontending 
node in iSAioha, transmits more frequently than a Type II node with the same transmission 
probability q 2 in Sd- Consequently, there are more eollisions in iSAioha eompared with those 
in 5/) and the serviee rate of Type I node in 5^:) eannot be smaller than the serviee rate 
of the first eontending node in iSAioha- Therefore, the stability eondition on Qi in iSAioha is 
suffieient for stability of Qi in Sd- 

(iii) In the stability analysis of Sd, we proved that the stability eondition on Qi is Ai < 
Qi (^1 “ T^)’ whieh is identieal to the stability eondition on Qi in iSAioha ETl - 

(iv) The suffieieney of the eondition Ai < gi for tho stability of Qi in Sd follows 

from its suffieieney for Qi in iSAioha- Therefore, the stability eonditions in (1^ are neeessary 
and suffieient eonditions and the region is the exaet stability region of the system Sd- 

■ 

E. The Relationship between Stability Regions of Sq and Sd 

Lemma 4. The stability region of the system Sd is an inner bound on the stability region of 
the system Sg, i-e.. Mg ^ M^. 

Proof: Aeeording to the assumption ¥[B > 0, Qi = OIQ 2 > 0] = 0 in the system Sd, a Type II 
node has a lower serviee rate eompared to a Type II node in the system So- Henee, the length of 
Q 2 in the system Sd is never smaller than its eounterpart in the system Sq, i.e., the system Sd 
dominates So from the perspeetive of Type II node. Consequently, the stability of Q 2 in Sd, is 
suffieient for the stability of Q 2 in Sg- Additionally, the assumption ]P[5>0,Qi = 0|Q2>0]=0, 
implies that Type II node is not transmitting in the time slots at whieh Qi is empty. Henee, 
the number of idle slots inereases, sinee at those time slots Type I node has no data paekets to 
transmit. Aeeordingly, the serviee rate of Type I node is not affeeted, and the stability eondition 
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on Qi is the same in 5/) and Sq- We eonelude that the stability region of is an inner bound 
on the stability region of iSg- ■ 

Finally, from Proposition 1 and Lemma 4 we arrive to Theorem 1. 


F. The Closure over all Transmission Probabilities 


In this subseetion, we prove Theorem 2. The elosure of the inner bound is defined 

by the union of all stability regions for a given (^1,^2), i-C-, C(r inner\ — u Mr^((gi,g 2 )). In 

(91,92) 

the serviee rate of Type II node is always lower than the arrival rate of Type I node, i.e., 
/i 2 < Ai- Also, from (fTSl) . we note that ^2 is inereasing in min{^^^, ^ 2 }, while /ii is not 
affeeted by 52 - Henee, in order to find the elosure we need to find q 2 that maximizes /i 2 - 

From (fTSl) . we observe that for maximizing /i 2 , the transmission probability of Type II node q 2 
should be greater than or equal to q^ = {qiPh)/{l + qiPh)- Also, inereasing q 2 beyond q^ does 
not affeet the serviee rate p 2 - 

Interestingly, we ean interpret q^ using Renewal reward theorem [[33l. Assume that the data 
queues are baeklogged and Type II node transmits whenever it reeeives an energy paeket, i.e., 
g 2 = 1- Let the expeeted reward R that Type II node obtains, be the transmission of one data 
paeket, i.e., E[i?] = 1. Also, the expeeted number of time slots, T, needed for the transmission 
of one data paeket is one slot for transmission, and {qiPh)~^ slots are needed for harvesting 
one energy paeket. Henee, the expeeted time needed for a transmission E[T] = 1 + {qiPh)~^- 
Using the renewal reward theorem, we find that the effeetive transmission rate of Type II node 
is given by g 2 ,e# = ||^ = Therefore, the previous expression represents the maximum 

possible transmission rate of Type II node, whieh is the minimum transmission probability q 2 
that maximizes yU. 2 . Now, the problem of finding the elosure C(§^), reduees to finding the elosure 
of ^ 2 )) over all gi, i.e., C(Rm„e.) = [j gj)), where 

9lG[0,l] 


j^inner 

G 


((91,^2)) 



Ai < gi(^l 


A2 \ 
1-qJ 


, A 2 < (1 - gi)pftAi 




(34) 


whieh represents the triangle OBD in Fig. where D = (gi, 0), and B = ■ 

Sinee, we know that the region eonsist of two line segments, the elosure C(Kinner) ean be 
found be taking the union of the elosures of the line segments OB, BD. First, we find the elosure 
of the line segment OB by solving x = and y = The solution represents the 
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trace of the point B for qi e [0,1], whieh is y = Phx{i. — Henee, the elosure of OB is 

given by _ ^. ^ 


~ I (-^ 1 ;-^2) A2 < 


AiPh 


(35) 


whieh is represented in Fig. [5] by Cq^ and Cq^. Note that C-^ is a eonvex region, sinee it 
is a hypograph of a eoneave funetion. Next, in order to find the elosure of BD, we solve 

max Qi (1 — ) • The solution gives us the elosure of the line segment extending OB to 

(/ie[o,i] V ^ _ 

the A 2 -axis, represented by and in Fig. [51 However, sinee we want the elosure of BD 
and not the extension, C-g^ is limited from the left by the trace of the point B. Thus, the elosure 
of BD is the region bounded from the left by and bounded by from the right. It is 
eharaeterized by 


^BD — 


.,A.) I 1. A.>A.p,(i-^-^) 


(36) 


2p{ ’ ' I-Aipfe. 

Note that the seeond eondition ean be found by solving the two equations of Cq^ and C\jj. 
Finally, the elosure is the union of Cqb and Cbd, 'whieh is represented by and Cgg, in 

Fig. [51 After some algebraie manipulations, we obtain (ITfil) . 

Remark: (ITfil) ean also be obtained by maximizing the serviee rate of Type II node /i 2 , subjeet 
to the stability eondition of Type I node Ai < pi, i.e. 


max (1 - qi)phXi s.t Ai < qi(l - )■ (37) 

9iG[o,i] V 1-gi/ 

Henee, the elosure of our system is equivalent to maximizing p 2 , beeause p 2 is upper bounded 
by pi. Thus, by maximizing p 2 , we implieitly maximize pi. 


VI. Model Extensions 

In this seetion, we diseuss two extensions to the previous stability analysis. First, we investigate 
the impaet of having a finite eapaeity battery on the stable throughput region. Seeond, we 
investigate the effeet of having a full-duplex RF energy harvesting node. 


A. Finite Capacity Battery 

In this subseetion, we investigate the impaet of having a finite eapaeity battery on the stable 
throughput region obtained in Theorem 1. Let M be the eapaeity of Type II node battery. Thus, 
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Fig. 5: The closure inner^ ~ ^OB Ph — 0-2. 


the battery evolution equation becomes 

Bt+i 

= min —/ig + if*, M} (38) 


Similar to the unlimited battery capacity case, the battery queue forms a decoupled discrete time 
Markov chain given that the data queues are backlogged. By analyzing the Markov chain we 
find the probability that the battery queue is non-empty. 

Lemma 5. For a half-duplex RF energy harvesting node with battery of size M, the probability 
that the battery is non-empty, C, given the data queues Qi and Q 2 are backlogged, is given by 


c 


p 



' 9iPh.(i-g2) 

92 



, 92 ^ 


<tiPh 
i+gip/i ’ 


92 


giPh 

^+qiph 


(39) 


where p = is the probability that the battery is non-empty in the infinite capacity 

battery case. 

Proof: Along the lines of Lemma 2. ■ 

Next, we apply the same procedure used in proving the stability region for the infinite battery 
capacity case. 

Corollary 1. An inner bound on the stable throughput region of the opportunistic RF energy 
harvesting slotted Aloha is the triangle OED, shown in Fig Assuming half-duplex energy 
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Fig. 6: The stable throughput region of opportunistie energy harvesting Aloha with finite 
eapaeity battery is eharaeterized by the triangle OED. 


harvesting and a battery of size M, the region is eharaeterized by the following equation 

K" = {(A. A.)| A. < (l - . A. < (40) 

Proof: Along the lines of Theorem 1. ■ 

The reduction in the stability region due to finite capacity battery is the triangle OEB shown 
in Fig. [6l 

B. Full-Duplex Energy Harvesting 

Now, we investigate the effect of having a full-duplex RF energy harvesting Type II node, i.e., 
the harvesting circuit is separated form the transmission circuit. Hence, a node can transmit and 
harvest simultaneously. Also, full-duplex energy transmission is advantageous because harvesting 
self-interference may introduce high energy yield. In the full-duplex RF energy harvesting 
paradigm, we have three types of harvesting opportunities. First, harvesting the transmissions of 
Type I node while Type II node is silent. Similar to the half-duplex case, we model this case by 
a Bernoulli process with mean Ph\{i}- Second, harvesting the self-interference of Type II node 
while Type I node is silent, which is modeled by a Bernoulli process with mean Ph\{2}- Third, 
harvesting both transmissions of Type I node and the self-interference of Type II node, which 
is modeled by a Bernoulli process with mean ph\{i,2}- 

In order to characterize the probabilities Ph\{2} and Ph\{i,2}, we use a similar approach to 
the one used in characterizing Ph\{i} in Section Hill We assume that the loopback interference 
coefficient c e [0,1] is known [|34]| . Also, we assume a Rayleigh fading channel between the 
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transmit antenna and the harvesting antenna, i.e., gt ~ exp(l). In case only self-interference is 
present, the received power at time slot t is equal to Pnit) = rjP 2 cgt. Hence, using the same 
approach as in the half-duplex case, we obtain 


1 

Vh\{2} = ^ 

^ VP2C 


(41) 


In case both transmissions of Type I node and the self-interference are present, the received 
power at time slot t is equal to Pnit) =viPihti^ + P) ^ + P 2 c gt). The probability Ph\{i, 2 } 
can be characterized in a similar fashioit. 

For a full-duplex harvesting node, the average energy harvesting process of the battery queue 
is given by 


K[H*'] — giPh\{i}^[Qi > 0] -f q2Ph\{2}^[Qi — 0, i? > 0, (^2 > 0] (42) 

+ {q2Ph\{2} + qiq2{ph\{i,2} - Phm - Ph|{i}))P[Qi > 0 ,5 > 0, Q2 > 0], 


where Ph\M^ is the harvesting probability given a set Ai of nodes are transmitting. The battery 
queue forms a decoupled Markov chain given that the data queues are backlogged. By analyzing 
the Markov chain we find the probability that the battery is non-empty, which is given by the 
following lemma. 

Lemma 6. For a full-duplex RF energy harvesting node with infinite capacity battery, the 
probability that the battery is non-empty \1/, given the data queues Qi and Q 2 are backlogged, 
is given by 


\1/ = min 


qiPh\{i} 


(43) 


, ^ 2(1 - qi{Ph\{i, 2 } - Ph\{i}) - (1 - qi)Ph\{ 2 }) ’ 

Proof: Along the lines of Lemma 2. ■ 

We notice that the probability of non-empty battery for the full-duplex case is higher than 
that of the half-duplex case, i.e, \1/ > The stable throughput region of the system is 

given by 

Corollary 2. An inner bound on the stable throughput region of the opportunistic RF energy 
harvesting slotted Aloha, under full-duplex energy harvesting mode and infinite capacity battery. 


^In order to characterize Ph\{i, 2 }. we need the distribution of the sum of independent gamma distributed random variables, all 
with integer shape parameters and different rate parameters, which is called the generalized integer gamma distribution (GIG) 
(35). 
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(a) q 2 > 


(b) 


91Ph|{l} 

l+91Ph|{l} 


< 92 < 5 '. 


(c) 92 < 


(llPh\{l} 
^+qiPh\{i} ■ 


Fig. 7: The stability regions of RF EH-Aloha under half/full-duplex modes vs. slotted Aloha 

with unlimited energy supply. 


is the region eharaeterized by 


Ig^inner 

G 


|(Ai, A2)| Ai < qi 


^2 \ 72(1 — 7 i)^Ai I 

1-gJ’ 71(1-72^) J 


Proof: Along the lines of Theorem 1. 


(44) 


VII. Numerical Analysis 
A. RF EH-Aloha System Sq Vs. Slotted Aloha 

In Fig. |71 we eompare the stable throughput region of the eonventional slotted Aloha with 
unlimited energy supply with our RF EH-Aloha system Sq- The stability regions are shown for 
qi = 0.4, Ph\{i} = 0.2, Ph\{ 2 } = 0.2, and Ph\{i, 2 } = 0.35. Also, we eonsider different values for q 2 
to eompare between full-duplex and half-duplex energy harvesting. We observe that the stability 
region of slotted Aloha is signifieantly redueed when RF energy harvesting is implemented, due 
to the energy limited sub-region. Also, for small q 2 , i.e., q 2 < observe that the 

stability regions of half-duplex and full-duplex EH-Aloha are identieal, beeause the serviee rate 
of Type II node is limited by the transmission probability q 2 . On the other hand, for large q 2 , 
i.e., q 2 > 4^, full-duplex RF energy harvesting expands the stability region, whieh agrees with 
intuition. From Fig. 13a), we notice that increasing q 2 beyond 4^ enhances the throughput of 
node 2 in the slotted Aloha system. However, the throughput of the Type II node in full-duplex 
EH-Aloha is limited by 4^. 
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Fig. 8: The stability regions 'Rd illustrated with the simulated sum average queue lenghts. 


B. 'Validating the Stability Conditions using Simulations 

In Fig. [8] and HI we simulate the energy harvesting Aloha system So for qi = 0.4, q 2 = 0.4, 
Ph = 0.6 and the queues are averaged over 10^ time slots. In partieular. Fig. [8] shows the sum of 
the average queue lengths E[Ql]+E[Q 2 ] versus the arrival rates (Ai, A 2 ). We observed that the 
system exhibits unstable behavior, shown by the inerease in the average queue lengths, as we 
eross the boundaries of the stability region. Fig. |9] shows the average serviee rate of Q 2 versus 
the arrival rates (Ai,A 2 ). It shows that the maximum serviee rate (lower left eomer points in 
the eontour plot), for a given A 2 , is aehieved on the boundary of the stability region. These 
observations support that in Proposition 1 is indeed the stability region of Sd- 

For the aforementioned parameters, we simulate the system So in order to verify the inner 
bound in Theorem 1. Similarly, the sum of the average queue lengths and the average serviee 
rate of Q 2 are shown in Fig. [TO] and dH respeetively. We notiee from Fig. (TOl that there exist 
rate pairs outside the left hand side of the stability region for whieh the queues exhibit a stable 
behavior. Additionally, Fig. [H] suggests that the maximum serviee rate (lower left eorner points 
in the eontour plot), for a given A 2 , is aehieved outside the stability region. These observations 
indieates that is an inner bound on the stability region of Sq, as proposed in Theorem 1. 

Finally, we show the utility of the stability region for the exaet system So whieh is 

deseribed in Section IlII-Al We simulate the exact behavior of the system for r] = 0.7, 7 = 0.2335, 
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Fig. 9: Simulation of the average service rate of Q2 in system Sd- 
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Fig. 10: The stability region illustrated with the simulated sum average queue lenghts. 


Pi = 1 and (1 + /")“^ = 0.5. Hence, 0 = 0.667 and from dH) we get ph = 0.6. The sum average 
queue lengths is shown in Fig. [I2l versus the arrival rates (Ai, A 2 ). We observe that the behavior 
of the average queue lengths in So is very similar to Sq- Henceforth, we conjecture that the 
stability region is also an inner bound for the stability region of the exact system Mq. 

To further support our analytical results. Fig. [13] and [Ml present sample paths for the evolution 
of the queues Qi and Q 2 in the systems So, So and Sq- Note that the sample path of the 
evolution of an unstable queue should show an increasing tendency such that the queue size 
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(a) 3D plot 


Fig. 11: Simulation of the average service rate of Q 2 in system Sq- 
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Fig. 12: The stability region compared with the simulated sum average queue lenghts of 

the exact system So- 


grow unboundedly as time increases. In Fig. [T3l we show the evolution of Q 2 for two arrival 
rate pairs Tq and Tj. We observe that Q 2 exhibits unstable behavior at Tq in So, while it exhibits 
stable behavior in So and Sq- This supports our claim that the stability condition on Q 2 in 
(flSl) is sufficient and necessary in Sd, while it is only sufficient in So and So- While, Fig. [I4] 
suggests that the stability condition on Qi in (fT5l) is sufficient and necessary in So, Sq and So- 
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Fig. 13: Sample paths for the evolution of Q 2 in the systems So, Sq and So- 



Fig. 14: Sample paths for the evolution of Qi in the systems Sd, Sq and Sq- 


VIII. Conclusion 

In this paper, we studied the effeets of opportunistie RF energy harvesting on the stability of 
a slotted Aloha system eonsisting of a Type I node, whieh has unlimited energy supply, and a 
Type II node, whieh is solely powered by an RF energy harvesting eireuit. We illustrated the 
intrieaey in analyzing the exaet behavior of sueh systems and proposed an equivalent system for 
whieh we were able to derive analytieal results. In partieular, we eharaeterized an inner bound 
on the stability region under the half-duplex and full-duplex energy harvesting paradigms, by 
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generalizing the stochastic dominance technique for RF EH-networks. We verified our analytical 
findings by simulating the exact and equivalent systems. The extension of our analysis to a 
random access network with multiple nodes can provide further insights to the development 
of efficient medium access protocols for networks with RF energy harvesting capabilities, and 
presents itself as a promising future research direction. 
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